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NOMENCLATURE 

aj  constant  defined  by  eq  3  where  i  =  0, 1 
A  constant 

Ai  constant  where  t  =  0, 1 

di  density  of  the  tth  constituent 

e  void  ratio 

/{  mass  flux  of  the  ith  constituent  relative  to  that  of  soil  minerals  where  i-1,2 

fiQ  mass  flux  of  water  in  the  unfrozen  part  of  the  soil 
h  amoimt  of  heave 

ko  thermal  conductivity  of  the  unfrozen  part  of  the  soil 

ki  thermal  conductivity  of  an  ice  layer 

Kj  empirical  function  defined  by  eq  8a  where  i  =  1, 2 

Kji  limiting  value  of  as  jc  approaches  n|  while  x  is  in  Rj,  t  =  1,2 

L  latent  heat  of  fusion  of  water,  334  Jgr^ 

m  location  of  the  top  end  of  the  specimen 

Mi  name  of  a  model  where  i  =  1, 2, 3 

n  boundary  in  where  the  pressiue  of  water  is  specified 

«;  boimdaiy  with  i  =  0, 1  where  hq  denotes  the  boundary  where  T  =  O^C  and  nj  the 

interface  between  an  ice  layer  and  a  frozen  fringe 
Pj  pressure  of  the  ith  constituent  where  /  =  1, 2 

Pio  value  of  Pj  at  tio 

^21  value  of  P2  at  ni 

r  rate  of  frost  heave 

Ro  imfrozen  part  of  the  soil 

R]  firozen  fringe 

R2  ice  layer 

R„  region  in  the  diagram  of  temperature  gradients  where  an  ice  layer  melts 

Rg  region  in  the  diagram  of  temperature  gradients  where  the  steady  growth  of  an 

ice  layer  occurs 

RJ  boundary  between  R,  and  Rj, 

RJ*  boimdary  between  R„  and  Rg 

Ru  region  in  the  diagram  of  temperature  gradients  where  the  steady  growth  of  an 

ice  layer  does  not  occur 
S  defined  by  eq  10 

So  value  of  S  when  o  =  0 


IV 


SPo  defined  by  eq  2 

t  time 

T  temperature 

T 1  temperature  at 

Ti*  temperature  at when  eq  1  holds  true 

Tf.  temperature  of  the  cold  bath 

Ty,  temperature  of  the  warm  bath 

X  spatial  coordinate 

Oq  absolute  value  of  the  temperature  gradient  at 

tt]  absolute  value  of  the  limiting  temperature  gradient  as  x  approaches  tii  while  x  is 
in  R2,  defined  by  eq  5 

Of  absolute  value  of  the  temperature  gradient  near  nj  in  R2 

ocu  absolute  value  of  the  temperature  gradient  near  Hq  in  Po 

y  constant,  1.12  MPa 

o  effective  pressure  defined  by  eq  1 

i  subscript  denotes  the  tth  constituent  of  the  mixture  consisting  of  unfrozen  water 
(t  =  1),  ice  (i  =  2)  and  soil  minerals  (t  =  3) 

*  superscript  used  to  indicate  the  value  of  any  variable  evaluated  when  a  point  (ui, 
Oq)  in  the  diagram  of  temperahire  gradients  is  on  Rt 
**  superscript  used  to  indicate  the  value  of  any  variable  evaluated  when  a  point  (ai, 
Oq)  in  the  diagram  of  temperature  gradients  is  on 
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Growth  Condition  of  an  Ice  Layer  in  Freezing  Soils 
Under  Applied  Loads 

I.  Experiment 

KAZUO  TAKEDA  AND  YOSHISUKE  NAKANO 

INTRODUCTION 

When  moist  fine-grained  soil  freezes,  it  often  is  accompanied  by  volume  expansion  caused  by  the 
appearance  of  layers  of  more  or  less  pure  segregated  ice  vsdthin  the  soil.  The  scientific  investigation 
of  segregated  ice  began  in  the  early  1900s,  but  the  imderstanding  gained  by  Taber  (1929,  1930), 
Beskow  (1935)  and  others,  although  useful,  was  largely  qualitative.  Since  then  a  significant  amoimt 
of  effort  has  been  made  to  gain  a  quantitative  imderstanding  on  ice  segregation. 

We  will  consider  the  one-directional  steady  growth  of  an  ice  layer.  Let  the  freezing  process 
advance  from  the  top  down  and  the  coordinate  x  be  positive  upwards  with  its  origin  fixed  at  some 
point  in  die  unfrozen  part  of  the  soU.  A  transitional  zone,  often  referred  to  as  the  frozen  fiinge,  exists 
between  the  frost  front  (0®C  isotherm)  and  the  growing  surface  of  an  ice  layer,  though  the  nature  of 
this  zone  has  not  been  well  understood.  A  freezing  soil  in  this  problem  may  be  considered  to  consist 
of  three  parts:  the  unfrozen  part  Rq,  the  frozen  fringe  and  the  ice  layer  R2,  as  shown  in  Figure  1. 

Radd  and  Oertle  (1973)  studied  the  relationship  between  the  pressure  P21  and  the  temperature  Tj 

of  an  ice  layer  at  by  using  a  water-permeable,  con-  _ 

stant-volume  cell  in  which  the  unfrozen  part  of  a  soil 
column  was  connected  with  a  water  reservoir  subject¬ 
ed  to  an  atmospheric  pressure,  Pj  =  0.1  (MPa).  They 
foimd  that  there  is  a  unique  temperature  T*  *  for  given 
P21  and  Pj  when  an  existing  ice  layer  neither  grows  nor 

melts  and  the  flux  of  water/,  vanishes,  and  that  T**is  - T  = 

given  as: 

o  =  P2i-Pi  =  -YTr,  if/i  =  0  (1) 

where  a  and  P21  are  often  referred  to  as  the  effective 
pressure  and  the  overburden  pressure,  respectively, 
and  Y  is  a  constant  with  the  value  of  1.12  MPa  ®C->. 

Equation  1  is  often  called  the  generalized  Clausius- 
Clapeyron  equation,  which  is  attributed  to  Edlefeen 
and  Anderson  (1943).  Studying  the  relationship  be¬ 
tween  the  effective  pressure  o  and  Ti*  by  using  a 
closed  constant-voliune  cell,  T2ikashi  et  al.  (1981)  con¬ 
firmed  the  validity  of  eq  1. 


Konrad  and  Morgenstem  (1980, 1^1, 1982)  empirically  found  that  the  rate  of  water  intake,/io,  at 
the  formation  of  the  final  ice  lens  is  proportional  to  foe  average  tonperature  gradieit  (T),  in  the 
finzen  fringe  when  a  soil  sample  freezes  under  difiermt  cold-side  step  temperatures  but  dw  same 
warm-side  temperature  under  atmospheric  conditions.  This  may  be  written  as: 


/io  =  -SPo(n.  (2) 

where  a  prime  denotes  differentiation  with  respect  to  x.  The  positive  proportionality  factor  SPg  was 
fotmd  to  be  constant  for  a  given  soil  and  was  termed  foe  segregation  potential.  They  also  found 
empirically  that  SPg  depends  on  the  pressure  of  water  Pjq  at  and  foe  applied  pressure  a. 

Studying  foe  relationship  between  fw  and  (T')«  experimentally,  Ishizald  and  Nishio  (1985) 
confirmed  foe  validity  of  eq  2  at  foe  formation  of  foe  final  ice  lens.  They  also  foimd  that /^g  is  a  linear 
function  of  foe  temperature  difference,  AT =Ti*  -T^duringfoestablegrowfoperiodoffoefinalice 
lens,  namely: 


AT  =  Og  +  fli/ig 

(3) 

AT=Tf*-Ti 

(4) 

where  <ig  and  are  constants.  Nakano  and  Takeda  (1991)  confirmed  empirically  foe  validity  of  eq  3 
when  0  =  0. 

We  studied  mathematically  and  experimetttally  foe  steady  growth  condition  of  an  ice  layer  under 
negligible  applied  pressure  in  three  earlier  reports  in  this  series  (Nakano  1990,  Takeda  and  Nakano 
1990,  Nakano  and  Takeda  1991).  We  have  shown  that  for  a  given  hydraulic  condition  foe  steady 
growth  condition  is  determined  by  foe  absolute  value  of  foe  temperature  gradient  Og  at  rtg  and  foe 
limiting  value  of  foe  temperature  gradient  in  R2  defined  as: 


ai  =  -  lim  T'(x).  (5) 

ztni^ 


The  steady  growth  condition  is  presented  in  foe  diagram  of  temperature  gradients  as  shown  in 
Figure  2.  The  steady  growth  occurs  in  foe  region  boimded  by  RS*  and  R,*.  The  boundary  RS*  is 


a  straight  line  given  as 

ao  =  (k,/fco)ai  (6) 

where  ki  and  are  foe  thermal  conductivities  of 
R2  and  Rg,  respectively.  We  will  refer  to  foe 
region  as  R„  where  Og  >  (ki/lfl)ai,/ig  <  0  and  an 
ice  layer  is  melting.  It  should  ^  mentioned  that  a  g 
an  existing  ice  layer  neither  grows  nor  melts,/ig 
vanishes  and  eq  1  holds  true  on  R**. 

The  boundary  R/  is  approximately  a  straight 
line  given  as: 

ao  =  Soai  ki  kg'  >  So  >  0  (7) 

where  Sg  is  a  constant  that  depends  on  only  foe 
properties  of  a  given  soil  if  ag>2.0  (®C  cm-*).  The 
steady  growth  of  an  ice  layer  does  not  occur  in 
foe  region  R^  where  <Xg  <  SgOj. 


Figure  2.  Temperature  gradients  and  Og. 
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WemtioduoedamodelMi(Nakano  1990)  ofafrozen  fringe  where  kemayexist  but  does  not  grow. 


and  the  mass  flux  of  water  fi(-  /iq)  is  given  as 

=  forxinRi 

dx  dx 

(8a) 

K2/K1  -4  y  as/j  ->  0 

(8b) 

lim  Pi(a^  =  P2(ri)  =  Pa 

(8c) 

xmki 

where  Ki  and  K2  are  properties  of  a  given  soil  that  generally  depend  on  T  and  the  compositian  of  die 
soil. 

In  the  three  earlier  reports  in  diis  series  we  showed  diat  is  consistent  with  eq  1  (N akano  1990), 
eq  3  and  7  under  o  =  0  (Nakano  and  Takeda  1991).  We  also  showed  (Nakano  and  Takeda  1991)  that 
the  boundary  R}  under  o  =  0  is  given  as: 

oo  =  ki{ko  +  LbKnV  ai  (9) 

where  L  is  the  latent  heat  of  fusion  of  water,  bis  a  function  of  the  thkkness  6  of  Ri  defined  by  eq  73c 
and  73e  in  Nakano  (1990),  and  X21  is  the  limiting  value  of  X2  ss  x  approadies  R]  while  x  is  in  R] .  An 
asterisk  denotes  that  R21  is  the  value  of  R21  when  a  point  (oq,ao)  is  on  R/  in  die  diagram  of  temperature 
gradients. 

This  is  the  first  of  a  two-part  presentation  on  the  effects  of  applied  pressure  (XI  die  steady  growth 
c(xidition  of  an  ice  layer.  The  objective  of  this  work  is  to  show  empirically  how  the  applied  pressure 
affects  the  steady  growdi  ccmdition  of  an  ice  layer.  In  the  following  work  we  will  show  diat  the  model 
Ml  is  cxmsistent  widi  experimental  data  presented  in  this  work.  Moreover,  we  wUl  show  diat  Mi  is 
consistent  with  the  reported  empirical  relatirxiships,  such  as  eq  2  and  3,  under  applied  pressure. 

EXPERIMENT 

The  same  test  apparatus  used  in  the  previous  study  (Takeda  and  Nakano  1990)  was  ccxnbined  with 
a  loading  device  that  imposed  a  specified  effective  pressure  o  cm  a  test  specimen  as  shown  in  Figure 
3a.  Since  the  weight  of  this  loading  device  is  balanced  by  the  counterweight,  only  a  specified  and 
constant  pressure  a  is  applied  on  the  moving  top  surfoce  of  the  specimen,  provided  that  the  fricticxi 
of  the  suspending  wire  for  the  counterweight  is  negligibly  smaU.  However,  die  fii(ri(xi  of  the  wire 
is  not  negligible  when  o  is  small.  Therefore,  a  weight  was  applied  direcdy  (xi  the  top  of  die  specimen 
when  a  was  less  dian  50  kPa.  Figure  3b  is  a  photograph  of  die  test  apparatus. 

A  schematic  drawing  of  the  test  cell  is  ^own  in  Figure  3c.  The  test  ceU  used  in  this  study  is 
essentially  the  same  as  diat  used  in  Takeda  and  Nakano  (1990)  except  for  minor  mcxlificaticxis 
described  below.  In  order  to  improve  the  accuracy  of  temperature  measurements,  cxie  more 
therm(Kouple  was  added  at  x  =  -1.0  cm  in  the  center  of  the  unfrozen  part  of  die  specimen  where  a 
one-dimensi(xial  coordinate  xis  intrcxluced  with  its  origin  at  the  level  of  die  water  reservoir,  as  shown 
in  Figure  3a.  Hence,  four  thermocouples  were  placed  at  x  =  0.0,  -0.5,  -1.0  and  -1.8  cm  in  the  center 
of  the  specimen  and  four  thermocouples  at  x  =  0.0, 05, 1.0  and  1.8  cm  along  the  inside  waU  of  the 
sample  holder.  In  order  to  suppress  the  heat  transport  through  die  loading  arm,  a  plastic  disk  was 
placed  (XI  die  top  of  an  aluminum  disk.  The  diameter  of  these  two  disks  was  slighdy  less  than  the 
inner  diameter  of  the  test  cell  so  diat  the  disks  smcxidily  slided  upward  during  foe  test 

In  Figure  3a  cxily  (xie  soil  specimen  is  shown.  However,  in  aU  freezing  tests  two  duplicated  soil 
specimens  were  tested  simultaneously  under  foe  same  ccxidition  to  ccxifirm  foe  reproducibility  of 
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а.  Schematic  of  test  apparatus. 

1.  Soil  spedmen. 

2.  Load. 

3.  Loadngaim. 

4.  Dial  gauge. 

5.  Offset  bar. 

б.  Upper  load-guide  bar. 

7.  Lower  load-guide  bar. 

8.  Water  reservoir. 

9.  Counterweight. 


c.  Schematic  of  test  cell. 

1 . Plastic  disk.  1 1  .Glass  beads. 


2.  Upper  cylinder. 

3.  Aluminum  disk. 

4.  Frozen  soil. 

5. Teflon  tape. 

6.  Thermocouple. 

7.  Sample  holder. 

8.  Unfrozen  soil. 

9.  Loading  arm. 
lO.O-rIng. 


12.  Lower  cylinder. 

13. Cold  bath. 

U.PosWon  of  the  top  surface  of  a 
soil  specimen  before  freezing. 

15.  Freezing  front. 

16.  Initial  water  head. 

17. Thermal  insulator. 

18.  Warm  bath. 

1 9.  Direction  of  water  flow. 


b.  Test  apparatus. 


Load 

♦ 


Figure  3.  Test  apparatus  and  test  cell. 


Figure  4.  Consolidation  characteristics  (e  vs. 
akPa)  qfKanto  loam. 


'"S  10  20  SO  100  200  500 

a.  Pressure  (kPa) 

test  results.  Kanto  loam  was  used  as  a  soil  specimen.  The  properties  of  this  soil  were  described  in  the 
previous  paper  (Takeda  and  Nakano  1990).  A  sample  soil  was  mixed  with  water  to  make  a  de-aired 
slurry.  The  slurry  was  placed  in  a  sample  holder  and  was  consolidated  in  steps  to  a  final  pressure  400 
kPa  in  10  days.  Hien  the  consolidated  specimen  was  kept  under  a  given  test  pressure  a  (<  390  kPa) 
for  more  than  48  hours  to  become  equilibrated. 

The  data  on  the  relationship  between  the  void  ratio  e  and  the  applied  pressure  o  in  the  rebound 
process  are  presented  in  Figure  4.  From  the  data  we  find  that  the  variation  of  e  in  the  whole  range  of 
<T  <  400  kPa  amounts  to  2%  at  most.  Although  the  variation  of  e  itself  is  negligibly  small,  some  of 
physical  properties  such  as  the  hydraulic  conductivity  may  be  affected  significantly  by  such  a  small 
variation.  We  will  further  discuss  the  effects  of  e  in  the  following  report. 

The  temperature  of  the  cold  bath  and  the  temperature  of  the  warm  bath  were  kept  constant 
at  specified  values  for  more  than  12  hours  so  that  the  temperature  profile  in  the  specimen  would  reach 
a  thermal  equilibrium.  The  specified  values  of  were  between  0°C  and  -2°C.  The  freezing  test  began 
with  the  ice  nucleation  induced  by  putting  a  very  small  piece  of  ice  on  the  top  surface  of  the  specimen. 

At  the  beginning  of  the  test  the  freezing  front  n^  rapidly  penetrated  into  the  specimen.  By  changing 
the  temperatures  and/or  T„  in  steps,  we  could  stop  the  freezing  front in  a  specified  region,  0 
<  X  <  1.0  cm.  We  let  an  ice  layer  grow  in  the  specified  region  for  two  reasons.  First,  one  of  four 
thermocouples  for  measuring  temperatures  in  the  unfrozen  part  of  the  specimen  was  located  at  z  = 
0.0  cm.  Secondly,  the  hydraulic  condition  in  our  tests  was  specified  by  the  distance  Sq  between  and 
n,  where  n  is  defined  as  the  interface  between  the  glass  beads  and  the  specimen  and  the  pressure  of 
water  was  kept  at  the  atmospheric  pressure.  Therefore,  it  was  important  to  conduct  all  tests  under 
nearly  the  same  hydraulic  conditions.  Since  n^  was  intentionally  positioned  in  the  region  specified 
above,  the  measurements  of  temperature  in  R2  were  made  by  the  thermocouples  along  the  inside  wall 
of  the  sample  holder.  Although  the  differences  between  the  measured  temperature  at  the  center  of 
the  specimen  and  that  along  the  inside  wall  at  x = 0.0  cm  were  found  to  be  quite  small,  the  temperature 
measurements  in  were  more  acciuate  than  those  in  Rz  because  the  temperature  at  the  center  was 
closer  to  the  average  temperature  over  the  cross  section  of  the  specimen. 

As  soon  as  the  freezing  front  stopped  and  the  thermal  field  in  the  specimen  attained  a  steady 
condition,  an  ice  layer  emerged  (initiation  of  ice  segregation).  From  the  measured  temperature  profile 
we  calculated  a  pair  of  variables  (a* ,  aS )  where  the  asterisks  denote  the  values  of  a  pair  (oCf,  oc„)  at  the 
initiation  of  ice  segregation.  As  we  will  show  below,  the  measured  temperature  profiles  in  Rq  (or  R^ 
were  nearly  linear  when  the  steady  growth  of  an  ice  layer  took  place.  The  value  of  ol^  (or  c%)  was 
calculated  for  the  most  part  from  the  measured  temperatures  at  the  two  points  nearest  to  Mq  (or  n^) 
in  Rq  (or  R^. 
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Table  1.  Teat  conditions. 


As  discussed  in  Nakano  and  Takeda  (1991), 
the  amount  of  heat  transported  by  convection 
is  much  less  than  that  transported  by  conduc¬ 
tion  in  our  experiments  and,  as  the  result,  the 
temperatureprofiles  in  l^(orR2)  become  near¬ 
ly  linear.  After  obtaining  the  first  pair,  we 
dianged  and/or  such  that  the  freezing 
front  moved  a  little  and  a  new  ice  layer  emerged 
at  a  new  site  slightly  below  the  previous  site. 

Repeating  the  same  procedure,  we  obtained 
several  pairs  of  (a*,  aS)  in  die  specified  region  described  above.  A  series  of  tests  were  conducted 
under  various  applied  presstues.  The  test  conditions  of  all  tests  are  summarized  in  Table  1 . 


a 

(kPa) 

No.cf 

runs 

Tc 

(“C) 

(•C) 

No.  ^ pairs 
(otf.aS) 

0.0 

14 

-6.1  to -20.0 

53  to  19.6 

14 

8.1 

4 

-3310-14.9 

3.0  to  15.1 

5 

16.2 

3 

-1.6  to -14.3 

1.6  to  14.0 

8 

48.7 

3 

-13  to -13.5 

1.6  to  13.0 

6 

97.5 

2 

-1.4  to -163 

3.0  to  283 

9 

195 

2 

-13  to -17.1 

2.0  to  323 

9 

390 

2 

-2.8  to -18.6 

5.0  to  353 

6 

TEST  RESULTS 

Several  variables  were  measured  dur¬ 
ing  the  test.  The  temperature  at  the  loca¬ 
tion  of  each  thermocouple  was  measured 
at  least  twice  every  hour,  while  the 
amount  of  heave  h  and  the  location  of  the 
freezing  front  were  measured  at  least 
once  every  hour.  The  amount  of  water 
taken  away  from  the  reservoir  was  mea¬ 
sured  at  suitable  intervals,  depending 
upon  the  amouint.  The  accuracy  of  these 
measurements  was  described  in  the  pre¬ 
vious  study  (Takeda  and  Nakano  1990). 

We  will  show  the  behavior  of  mea¬ 
sured  variables  with  time  by  using  the 
set  of  data  obtained  under  the  applied 
pressure  of  195  kPa  as  an  example.  The 
records  of  Tj  and  are  shown  in  Figure 

5a.  Since  and  Ty,  were  changed  step¬ 
wise,  the  records  look  like  a  step  func¬ 
tion.  Despite  the  stepwise  change  of  T^ 
and  Ty,,  all  other  measured  variables  de¬ 
scribed  above  changed  more  or  less  con¬ 
tinuously  with  respect  to  time.  The 
amount  of  heave  h  and  the  position  of 
are  plotted  vs.  time  in  Figvue  5b. 

The  rate  of  heave  r  was  calculated 
from  the  measured  h{t)  as  the  average 
over  a  finite  time  interval,  while  the  rate 
of  water  intake /jq  was  calculated  from 
the  measured  amount  of  water  in  the 
reservoir  as  the  average  over  another 
finite  time  interval.  The  calculated  val¬ 
ues  of  r  and /lo  are  plotted  vs.  time  in  Fig¬ 
ure  5c  where  each  horizontal  line  is  the 
average  over  the  time  interval  corre- 


a.  TemperaturesTc°Camn^°CvsAhoursxvitha=  195kPa. 


b.  Measured  h  cm  and  nj  cm  os.  t  hours  with  a  =  195  kPa. 
Figure  5.  Measured  variables  vs.  time  for  a  =  195  kPa. 
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spending  to  the  line  length.  Similarly,  the  calculated  average  values  of  Of  and  tx^  from  ti\e  measured 
temperature  profiles  over  certain  variable  intervals  are  plotted  vs.  t  in  Figure  5d. 

The  measured  temperature  profiles  in  the  specimen  at 73.90, 96.45, 104.43, 1 19.85  and  144.75 hours 
are  presented  in  Figure  5e,  where  dark  circles  are  data  points  and  straight  lines  are  drawn  between 
two  neighboring  points.  It  is  easy  to  see  from  Figure  5e  that  each  solid  line  in  the  frozen  or  the 
unfrozen  part  of  the  specimen  is  nearly  linear.  This  implies  that  the  thermal  field  was  nearly  steady 
at  each  of  those  points  of  time  when  the  initiation  of  ice  segregation  occurred  and  a  pair  of  (a*,  aS) 
were  obtained. 

At  the  beginning  of  this  test,  and  T„  were  kept  comtant  at  -0.9°C  and  2.0°C,  respectively.  The 
freezing  firont  Wj  moved  downward  rapidly.  At  48.3  hours  was  decreased  to  -1.3°C  so  that  the 
penetrating  speed  of  would  be  slowed  while  was  kept  coirstant.  At  about  55.0  hours  the  speed 
of  «i  increased  again  and  the  rate  of  heave  r  began  increasing  to  attain  a  maximiun  at  61 .0  hours.  After 
61.0  hours  r  decreased  and  gradually  approached  x = 0.85  cm.  At  71 .4  hours  a  segregated  ice  layer 
without  any  visible  soil  particles  emerged  as  soon  as  Wj  ceased  to  move.  At  fiiis  time  r  became  nearly 
constant  and  the  values  of  and  remained  constant.  Under  such  a  steady  condition,  the  ice  layer 
kept  growing  while  Mj  remained  still  until  73.9  hours.  We  foimd  these  constant  values  of  Of  and  a,, 
to  be  a*  and  aS/  respectively.  In  other  words,  a  constant  pair  (a*,  aS)  was  on  R**. 

After  obtaining  the  first  pair  (a*,  a5),  we  changed  and  in  such  a  manner  that  nj  moved 
downward  slightly  and  a  new  ice  layer  emerged  at  a  new  location  slightly  below  the  previous  ice 
layer.  Repeating  the  same  procedure,  we  obtained  several  pairs  of  (a*,  aS).  As  shown  in  Figure  5d, 
it  took  a  few  hours  for  Of  and  to  stabilize  after  and/ or  were  changed  stepwise.  It  is  clear  from 
Figures  5c  and  5d  tha^  il'.e  values  of  r  and /jq  were  sensitively  affected  by  minor  variations  of  and 
Oy.  When  r  is  quite  small,  accurate  measurement  of  the  rate  of  water  intake  over  a  short  time  interval 
was  difficult.  Therefore,  the  mass  flux  of  water  was  usually  calculated  based  on  r  under  such  a  case. 

From  two  tests  under  the  condition  of  o  =  195  kPa,  we  obtained  nine  pairs  of  (a*,  aS)  that  are 
plotted  in  the  diagram  of  temperature  gradients  in  Figure  5f.  As  in  the  previous  study  (Takeda  and 
Nakano  1990),  we  will  seek  a  linear  approximation  to  the  relation  between  a*  and  aS  given  as: 

aS  =  S(a)  a*  (10) 

where  the  positive  proportionality  factor  S  is  assumed  to  depend  on  only  the  applied  pressure  a. 
Using  the  method  of  interval  estimate  (Fisz  1965),  we  have  found  from  the  data  set  that  the  value  of 
S(195)  is  1.856.  The  line  R*  (195)  passing  through  the  data  points  in  Figure  5f  is  given  by  eq  10  with 
a  =  195  kPa,  which  is  a  linear  approximation  to  the  boundary  R*  that  divides  the  steady  growth 
region  and  the  region  R^  where  the  steady  growth  of  an  ice  layer  does  not  occur. 

The  broken  line  R  J  •  (195)  in  Figure  5f  is  given  as: 

Ou  =  [kiicykolOf ,  a  =  195  kPa. 

We  found  empirically  that  the  thermal  conductivity  ko  of 
the  vmfrozen  part  of  a  specimen  is  hardly  affected  by  the 
applied  pressure  <T.  However,  a  affects  the  thermal  con¬ 
ductivity  fci  of  an  ice  layer  to  a  degree  that  is  minor  but  not 
negligible.  The  measured  average  densities  ^2  and  con¬ 
ductivities  ki  of  ice  layers  imder  various  a  values  are  pre¬ 
sented  in  Table  2.  It  is  easy  to  see  that  the  values  of  di 
when  <i  >  0  are  greater  than  when  o  =  0.  The  applied 
pressure  evidently  reduces  the  amount  of  space  occu¬ 
pied  by  air  in  ice  layers  and,  as  a  result,  the  values  of  k^ 
when  o  >  0  are  also  greater  than  when  a  =  0  as  shown  in 
Table  2. 


(11) 


Table  2.  Summary  of  test  results. 


a 

(kPa) 

S(a) 

dz 

(g  cm-3) 

k, 

J  (cm  s  °C)-' 

■  O 

0.0 

0.871 

0.860 

1.59E-2+ 

1.73 

8.1 

1.078 

0.898 

2.13E-2 

2.31 

16.2 

1.554 

0.895 

2.12E-2 

2.30 

48.7 

1.618 

0.907 

2.19E-2 

2.38 

97.5 

1.817 

0.906 

2.18E-2 

2.37 

195 

1.856 

0.907 

2.19E-2 

2.38 

390 

1.985 

0.904 

2.17E-2 

2.36 

’  *„  =  9.20xl0-U/(cm  s°C). 
tE-N=10-N 
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a,  a, 

a.a  =  8.12kPa.  b.a=16.2kPa.  c.a  =  48.7kPa. 


d.<J  =  97.5  kPa.  e.a  =  390  kPa. 


Figure  6.  Steady  growth  region  Rg. 

The  values  of  (fci/fco)  with  a  =  195  kPa  and  a  =  0.0  kPa  are  2.38  and  1.73,  respectively,  as  given  in 
Table  2.  Because  of  tins  difference  the  broken  line  RS*  (0)  given  by  eq  11  with  <j  =  0.0  lies  below  the 
line  Ra  *  (195)  as  shown  in  Figure  5g.  The  steady  growtii  region  i?s(195)  under  a = 195  kPa  is  tihe  region 
bounded  by  Rt  (195)  and  R*  *  (195)  while  that  under  a = 0.0  kPa  is  the  region  bounded  by  Rt  (0)  and 
Rt*  (0)  in  Figure  5g.  It  is  clear  from  this  figure  that  the  applied  pressme  of  195  kPa  significantly 
reduced  the  area  of  the  steady  growth  region. 

Using  the  same  test  procedure  as  described  above  for  the  case  of  a = 195  kPa,  we  conducted  a  series 
of  tests  tmder  various  applied  pressiues  to  determine  the  region  of  steady  growth  R,(a)  for  a  given 
o.  The  results  of  these  tests  are  presented  in  Figure  6a  (o  =  8.12  kPa),  6b  (16.2  kPa),  6c  (48.7  kPa),  6d 
(97.5  kPa)  and  6e  (390  kPa);  It  is  clear  from  Figures  5g  and  6  that  the  region  of  steady  growth  Rg(G) 
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touls  to  decrease  with  increasing  o.  The  values  of  S(o) 
givot  in  Table  2  are  plotted  vs.  a  in  Figure  7.  It  is  found 
from  Figure  7  that  the  value  of  S  increases  sharply  wid\ 
increasinga  when  ais  small  and  that  the  rateof  increase 
slows  down  as  o  becomes  greater.  A  curve  in  Figure  7 
is  one  example  of  the  approximate  presentation  of  data 
points  given  as 

S(o)  =  0.280  o0^«  + 0.856  (12) 

where  o  is  in  kilopascals. 

We  have  found  empirically  (Takeda  and  Nakano 
1990)  that  there  is  the  upper  boimd  Aj  (or  A^)  of  Of  (or  oj 
beyond  which  an  ice  layer  with  cavities  grows  when  o 
is  negligibly  small  and  that  the  value  of  Aj  for  Kanto 
loam  is  8.8°C  cm-^  The  reason  for  such  behavior  was 
explained  (Takeda  and  Nakano 1990,  Nakano  and  Take¬ 
da  1991)  as  follows.  The  pressure  water  at  ng  decreases  with  the  iiKreasing  value  of  otf  (or/jg) 

imtil  Pio  attains  a  value  that  corresponds  to  the  air  entry  value  of  Kanto  loam  when  an  ice  layer  wiA 
cavities  begins  growing.  The  growth  of  an  ice  layer  with  cavities  was  not  observed  in  all  the  tests 
under  various  applied  pressures  presented  in  this  work.  Consequently,  we  were  unable  to  confirm 
the  upper  bound  of  oq  in  this  work. 

CONCLUSIONS 

A  series  of  freezing  tests  on  Kanto  loam  were  conducted  under  various  overburden  pressures  to 
find  the  steady  growth  condition  of  a  segregated  ice  layer.  The  steady  growth  condition  was  foimd 
to  be  determined  by  the  absolute  value  of  the  temperature  gradient  of  fi\e  ui\fit>zen  part  of  fite  soil 
a„  near  the  0°C  isotherm  and  fiiat  of  the  frozen  part  of  the  soil  o({  near  the  warmest  end  of  an  ice  layer 
under  given  hydraulic  condition  and  applied  effective  pressure  as  follows: 

OusAo,,  lci/fco>  A  >S(o)  (13) 

where  S(a)  is  an  increasing  function  of  o. 

This  is  the  first  of  a  two-part  presentation  on  the  efiects  of  applied  pressure  on  the  steady  growfii 
condition  of  an  ice  layer.  In  the  following  paper  we  will  show  Aat  the  model  is  consistent  wifii 
experimental  data  presented  in  this  work. 
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